UNCLASSIFIED 


AD 

AD-E403  737 


Technical  Report  ARWSE-TR-1 5036 


REDUCING  PEAK  POWER  IN  AUTOMATED  WEAPON  LAYING 


Joshua  Stapp 
Matthew  Tomik 


February  2016 


U.S.  ARMY  ARMAMENT  RESEARCH,  DEVELOPMENT  AND 
ENGINEERING  CENTER 

Weapons  and  Software  Engineering  Center 

Picatinny  Arsenal,  New  Jersey 


Approved  for  public  release;  distribution  is  unlimited. 


UNCLASSIFIED 


UNCLASSIFIED 


The  views,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author(s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation. 

The  citation  in  this  report  of  the  names  of  commercial  firms  or  commercially 
available  products  or  services  does  not  constitute  official  endorsement  by  or 
approval  of  the  U.S.  Government. 

Destroy  this  report  when  no  longer  needed  by  any  method  that  will  prevent 
disclosure  of  its  contents  or  reconstruction  of  the  document.  Do  not  return 
to  the  originator. 


UNCLASSIFIED 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-01-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  the  burden  to  Department  of  Defense,  Washington  Headquarters  Services  Directorate  for  Information  Operations  and  Reports  (0704-0188), 
1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any 
penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 

February  2016 


2.  REPORT  TYPE 

Final 


3.  DATES  COVERED  {From -To) 


4.  TITLE  AND  SUBTITLE 


REDUCING  PEAK  POWER  IN  AUTOMATED  WEAPON  LAYING 


5a.  CONTRACT  NUMBER 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


6.  AUTHORS 

Joshua  Stapp  and  Matthew  Tomik 


5d.  PROJECT  NUMBER 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 
U.S.  Army  ARDEC,  WSEC 
Weapons  Systems  &  Technology  Directorate 
(RDAR-WSW-I) 

Picatinny  Arsenal,  NJ  07806-5000 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

U.S.  Army  ARDEC,  ESIC 

Knowledge  &  Process  Management  (RDAR-EIK) 

Picatinny  Arsenal,  NJ  07806-5000 


10.  SPONSOR/MONITOR’S  ACRONYM(S) 


11.  SPONSOR/MONITOR’S  REPORT 
NUMBER(S) 

Technical  Report  ARWSE-TR-1 5036 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 


In  artillery,  the  process  of  aiming  a  weapon  is  referred  to  as  gun  laying.  This  report  describes  a  method 
to  calculate  motion  profiles  that  reach  a  given  lay  within  the  least  amount  of  time  while  reducing  the  amount  of 
peak  power  required  and,  therefore,  minimizing  the  forces  caused  by  acceleration. 


15.  SUBJECT  TERMS 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 
ABSTRACT 

18.  NUMBER 
OF 

19a.  NAME  OF  RESPONSIBLE  PERSON 
Joshua  Stapp 

a.  REPORT 

U 

b.  ABSTRACT 

u 

c.  THIS  PAGE 

u 

SAR 

PAGES 

16 

19b.  TELEPHONE  NUMBER  (Include  area 

code)  (973)  724-3051 

Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


UNCLASSIFIED 


UNCLASSIFIED 


CONTENTS 

Page 

Introduction  1 

Methods,  Assumptions,  and  Procedures  2 

Conventions  and  Variable  Definitions  2 

Fundamental  Equations  3 

Calculating  Time  Periods  3 

Calculating  New  Acceleration  Values  5 

Results  and  Discussions  7 

Conclusions  10 

Distribution  List  1 1 

FIGURES 

1  Trapezoidal  motion  profile  1 

2  Triangular  motion  profile  2 

3  Trapezoidal  time  adjustment  5 

4  Invalid  trapezoidal  time  adjustment  6 

5  Simulation  of  15  deg  motion  8 

6  Simulation  of  45  deg  motion  8 

7  Simulation  of  45  deg  motion  (extended)  9 


Approved  for  public  release;  distribution  is  unlimited. 

UNCLASSIFIED 


UNCLASSIFIED 


INTRODUCTION 

In  artillery,  the  process  of  aiming  a  weapon  is  referred  to  as  gun  laying.  Gun  laying  involves  a 
set  of  actions  to  align  the  axis  of  a  gun  barrel  so  that  it  points  in  the  required  direction.  This  alignment 
is  specified  relative  to  the  horizontal  and  vertical  planes.  A  gun  is  traversed  in  the  horizontal  plane 
and  elevated  in  the  vertical  plane  to  range  it  to  the  target.  The  traverse  and  elevation  values  make  up 
the  aiming  portion  of  the  firing  solution.  In  an  automated  weapon  system,  the  rotation  of  these  two 
axes  is  performed  using  actuators. 

A  typical  motion  profile  follows  the  velocity  versus  time  graph  shown  in  figure  1 .  An  object  will 
accelerate  at  a  given  rate  up  to  a  maximum  velocity  and,  after  a  period  of  time,  will  decelerate  at  a 
given  rate  until  stopping  at  the  target  position.  In  this  figure,  displacement  (or  distance)  is 
represented  as  the  area  under  the  curve,  and  acceleration  is  represented  as  the  slope  of  the  line.  In 
some  cases,  the  required  distance  may  be  short  enough  that  the  maximum  velocity  is  never 
reached.  This  condition  is  shown  in  figure  2.  In  the  case  of  laying  a  gun,  there  would  be  a  separate 
profile  for  both  the  traverse  and  elevation  axis.  Distance,  velocity,  and  acceleration  are  replaced  by 
angle,  angular  velocity,  and  angular  acceleration. 


Figure  1 

Trapezoidal  motion  profile 
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Typical  actuator  controllers  are  configured  by  setting  acceleration  and  maximum  velocity 
values.  When  commanded  to  move  a  given  displacement,  the  controller  executes  a  motion  profile 
using  those  stored  values.  The  values  used  are  determined  based  on  a  number  of  factors  including 
available  power,  maximum  motor  speed,  maximum  safe  slewing  speeds,  peak  available  motor 
torque,  and  desired  maximum  laying  time.  These  precalculated  values  guarantee  that  the  system 
remains  within  safe  operating  limits  while  executing  within  the  required  time  for  a  worst  case 
targeting  situation.  If  these  values  are  not  optimized,  most  laying  motions  will  result  in  one  axis 
arriving  at  its  commanded  angle  before  the  other.  Since  the  weapon  is  only  layed  once  both  axes  are 
in  place,  there  is  an  unnecessary  amount  of  power  spent  accelerating  and  decelerating  the  closer 
axis.  By  reducing  the  acceleration  and  deceleration  rates  of  the  closer  axis  to  follow  a  motion  profile 
that  arrives  at  its  target  at  the  same  time  as  the  further  axis,  one  can  reduce  peak  power  and  the 
forces  imparted  on  the  system. 


METHODS,  ASSUMPTIONS,  AND  PROCEDURES 


Conventions  and  Variable  Definitions 

Before  describing  the  formulas  to  solve  the  aforementioned  problems,  it  is  necessary  to 
define  conventions  and  variables.  Table  1  defines  the  variables  used  in  the  subsequent  calculations. 
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Table  1 

Variable  definitions 


Variable 

Definition 

e 

Displacement  angle 

01 

Displacement  while  accelerating  to  peak  velocity 

02 

Displacement  while  accelerating  from  peak  velocity  to  0 

0slew 

Displacement  while  slewing  at  a  constant  velocity 

CUmax 

Maximum  angular  velocity  setpoint 

COP 

Peak  angular  velocity  of  triangular  motion  profile 

01 

Angular  acceleration  used  to  increase  velocity 

02 

Angular  acceleration  used  to  decrease  velocity  (decelerate) 

ti 

Time  spent  accelerating  to  peak  velocity 

t2 

Time  spent  accelerating  from  peak  velocity  to  0 

tslew 

Time  spent  slewing  at  a  constant  velocity 

ttarget 

Target  amount  of  time  to  complete  motion  in 

Fundamental  Equations 

Kinematics  equations  describe  the  motion  of  objects  in  space.  The  following  equations  will  be 
used  throughout  this  report  and  are  valid  for  constant  angular  acceleration  in  one  dimension. 


9  =  90  +  co0t  +  I/2  °T2 

(1) 

co  =  coq  +  at 

(2) 

co2  —  co02  +  2a(9  —  90 ) 

(3) 

When  working  with  velocity  versus  time  graphs,  it  is  often  simpler  to  solve  for  zero  by  thinking 
in  terms  of  the  area  under  the  curve.  For  periods  of  constant  velocity,  this  is  a  rectangle  whose  area 
is  equal  to 

9  —  cot  (4) 

For  a  period  of  constant  acceleration,  there  are  three  equivalent  equations  that  can  be  used 
depending  on  the  variables  of  interest.  These  equations  treat  the  period  of  deceleration  as  if  they 
were  accelerating  (hence  the  absolute  value  for  acceleration).  The  assumption  is  that  the  period  of 
acceleration  will  always  either  begin  or  end  at  a  velocity  of  zero. 

0  =  l/2  cot  =  a)2/2|a|  =  V2  l«lt2  (5) 


Calculating  Time  Periods 

The  first  step  in  calculating  time  periods  is  determining  the  shape  of  the  profile.  Given  0,  ow 
ai  and  02  it  can  be  determined  whether  a  trapezoid  or  triangle  shaped  motion  profile  will  be  used. 
This  is  done  by  first  assuming  a  triangular  profile  and  calculating  the  peak  velocity  (cop)  attained 
using  only  the  ai  and  02.  The  angle  to  travel  can  be  broken  into  two  parts.  The  first  occurs  during  the 
acceleration  period  and  the  second  during  the  deceleration  period.  For  the  purpose  of  finding  the 
angle,  the  value  of  02  can  be  treated  as  positive  (area  under  the  curve  of  a  triangle). 
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9  =  9±  +  92  =  1/2aitZ  +  1/2\a2\t2  (6) 

The  velocity  starts  at  zero  and  rises  to  a  peak  value  during  ti  (eq.  7).  During  h  it  starts  at  the 
same  velocity  and  ends  at  zero  (eq.  8).  These  two  formulas  are  set  equal  to  each  other  in  equation  9 
and  then  solved  for  t2  in  equation  10.  Note  that  the  negative  value  of  02  results  in  a  positive  value  for 
t2. 


co  —  co0  +  a xtx  —  0  +  a-tt-L  —  a xtx  (7) 

U final  =  0  =  coo  +  a2t2  =>  CO  =  -a2t2  (8) 


<Xih  =  ~a2t2 


(9) 


^2  — 


«2 


(10) 


The  result  of  plugging  the  equation  for  t2  into  equation  6,  solving  for  ti,  and  simplifying  is 


h 


(11) 


Multiplying  this  value  by  the  absolute  acceleration  rate  of  cm  provides  the  peak  angular 
velocity  for  a  triangular  motion  profile. 


00  p  — 


(12) 


If  (Op  is  less  than  comax  ,  the  motion  profile  is  triangular.  Otherwise,  the  motion  profile  is 
trapezoidal. 

For  a  triangular  profile,  all  the  work  has  already  been  completed  with  ti  and  t2  represented 
by  equations  10  and  1 1  respectively.  Since  cop  is  never  reached,  tsiew  is  set  to  zero. 

If  the  shape  of  the  motion  profile  is  trapezoidal,  the  times  can  be  solved  for  using  equations  2 

and  5. 


ti  = 


ty  — 


a2 


t slew 


<9-6>i-6/2  _  e~ 


/ 2  _ 


ti  +  t2 


(13) 

(14) 

(15) 
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Calculating  New  Acceleration  Values 

Once  the  time  periods  have  been  calculated,  it  can  be  determined  which  axis  takes  longer  to 
travel  using  the  given  acceleration  and  maximum  speed  values.  The  faster  axis  can  then  be  slowed 
down  to  match  the  travel  time  of  the  slower  axis.  The  new  target  travel  time  will  be  referred  to  as 
target.  The  first  step  in  calculating  the  values  for  the  new  motion  profile  is  to  determine  if  the  solution 
profile  will  take  the  form  of  a  trapezoid  or  a  triangle.  To  do  so,  the  change  in  time  between  target  and 
tatai  is  calculated. 


At  -  ^target  t total  ("16) 

If  At  is  less  than  tsiew,  the  motion  profile  of  the  solution  is  trapezoidal.  The  rationale  behind  this 
is  demonstrated  in  figures  3  and  4.  As  ai  and  02  are  decreased  (the  slopes  are  flattened),  an  equal 
amount  of  angular  displacement  (area)  is  both  added  and  removed,  keeping  the  total  displacement 
the  same.  However,  at  some  point,  the  acceleration  removal  results  in  a  total  erosion  of  the  available 
tsiew  time  therefore  requiring  a  triangular  profile  solution.  Obviously,  if  the  initial  profile  shape  is 
triangular,  tsiew  is  zero  and  the  solution  will  always  be  triangular  as  well. 


CD 


New  accelerations 


Lost  0  =  Added  0 

At2  Tsiew  ^  At1  +  At2 


Figure  3 

Trapezoidal  time  adjustment 
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Figure  4 

Invalid  trapezoidal  time  adjustment 

Once  the  shape  of  the  motion  profile  is  determined,  use  the  appropriate  set  of  calculations  to 
determine  the  new  values  of  ai  and  02. 

For  both  profiles,  it  is  assumed  that  the  ratio  between  acceleration  and  deceleration  will  be 
maintained,  resulting  in  equation  17.  Therefore,  written  as 


CTi  (Xi new 

a2 new 


(17) 


Triangular  Profile 

For  a  triangular  profile,  there  is  a  known  displacement  angle  0  and  a  given  target  time 
target.  New  acceleration  values  are  being  calculated,  which  can  be  called  amew  and  a2new.  The  tmew 
and  t2new  can  be  defined  as  the  new  times  spent  accelerating  and  decelerating.  The  sum  of  these  is 
the  total  target  time. 


t target  ~  new  T  ^2 new 


(18) 


Using  equations  10,  17,  and  18,  tmew  in  terms  of  ttargetcan  be  solved  as 

a2 

hnew  ~  ^-target 

Substituting  equation  19  into  equation  2  results  in 

^2 


(i)p  —  0Cinew 


(a2  -  aj ttarget 


(19) 


(20) 
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Finally,  equation  20  is  plugged  into  equation  5  and  solved  for  amew. 

_  20(a2  -  ax) 

“ln™  “  «2(tto 


(21) 


Once  ainew  is  determined,  equation  17  can  be  solved  for  Q2new. 


_  ^2^1  new 
^2 new  ~  ~ 

ai 


(22) 


Trapezoidal  Profile 

In  order  to  calculate  the  trapezoidal  profile,  the  change  in  time  for  the  acceleration 
and  deceleration  periods  needs  to  be  calculated.  These  two  changes  in  time  sum  to  the  total  change 
in  time  At  as  calculated  in  equation  16. 


At  =  Att  +  A  t2  (23) 

Then,  Abusing  equations  10,  17,  and  23  can  be  solved 


A  tx 


«2 

(a2-ai) 


At 


(24) 


For  the  trapezoidal  solution,  what  is  depicted  in  figure  3  is  calculated.  To  do  so,  equation  2  is 
written  using  Ati  as  calculated  in  equation  25. 


^ max  ~  ^lnewC^l  "F  2At^) 


(25) 


The  result  for  solving  for  a  Inew  is 


_  ^ max 

ainew  ~  (ti  +  2Ati) 


(26) 


With  amew  solved,  a2new  can  be  solved  using  equation  22  from  the  triangular  profile  solution. 


RESULTS  AND  DISCUSSIONS 

The  calculations  presented  in  this  report  were  entered  into  a  spreadsheet  and  various 
scenarios  were  simulated.  The  calculated  motion  profile  values  were  then  fed  into  models  for  the 
elevation  actuator  of  an  existing  weapon  system.  The  peak  power  values  of  the  modified  profiles 
were  compared  against  the  original  profiles.  For  the  original  profile,  acceleration  and  deceleration 
rates  of  25  deg/s2  were  used  with  a  maximum  angular  velocity  of  20  deg/s.  These  values  represent 
the  elevation  profile  limits  for  the  weapon  system. 

In  figure  5,  a  displacement  of  15  deg  was  performed.  The  original  profile  completed  this 
motion  in  ~1 .5  sec  using  a  triangular  motion  profile.  When  the  target  time  was  extended  to  4  sec,  the 
acceleration  and  deceleration  rates  were  decreased  from  25  to  3.75  deg/s2  with  the  profile  remaining 
triangular. 
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In  figure  6,  a  displacement  of  45  deg  was  performed.  The  original  profile  took  ~3  sec  to 
complete.  When  the  target  time  was  extended  to  4  sec,  the  rates  were  decreased  to  1 1 .25  deg/s2. 
The  profile  began  as,  and  remained,  trapezoidal  in  shape. 


Velocity  vs  Time 


Figure  6 

Simulation  of  45  deg  motion 
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In  figure  7,  the  same  displacement  of  45  deg  was  performed,  but  the  target  time  was 
extended  to  8  sec.  The  new  profile  decreased  its  rates  to  2.8  deg/s2,  and  the  resulting  profile  was 
transformed  from  trapezoidal  to  triangular. 


Figure  7 

Simulation  of  45  deg  motion  (extended) 

Simulations  of  an  existing  system  were  performed  and  the  following  scenario  was  examined. 
The  gun  is  pointing  at  0  deg  azimuth  and  15  deg  elevation.  It  gets  a  command  to  lay  at  180  deg 
azimuth  and  30  deg  elevation.  The  travel  times  are  calculated  and  compared.  The  traverse  axis  will 
take  8  sec  to  complete  its  motion  using  acceleration  and  deceleration  rates  of  20  deg/s2  and  a 
maximum  velocity  28  deg/s.  The  elevation  axis,  using  the  same  values  as  mentioned  previously  (25 
deg/s2  and  20  deg/s),  would  complete  its  motion  in  1 .55  sec.  The  traverse  axis  takes  longer  to 
complete  and,  therefore,  becomes  the  new  target  time  for  the  elevation  actuator.  Application  of  the 
algorithms  described  in  this  report  resulted  in  a  reduction  of  acceleration  from  the  original  25  to 
0.9375  deg/s2.  The  peak  power  used  by  the  elevation  axis  was  reduced  from  743  to  128  W.  The 
azimuth  power  remains  the  same  as  its  original  profile  is  maintained. 

Calculating  the  overall  benefit  to  a  particular  system  is  difficult.  The  peak  power  is  reduced  as 
a  function  of  the  difference  in  travel  times  between  the  axes.  If  they  take  the  same  amount  of  time, 
no  savings  are  realized.  From  a  mechanical  perspective,  there  is  no  energy  savings,  as  the  same 
amount  of  mechanical  work  is  being  performed  regardless  of  the  time  taken  to  perform  it.  When  the 
electrical  components  are  taken  into  account,  this  is  no  longer  the  case.  By  lowering  the 
acceleration,  the  torque  required  to  perform  the  acceleration  is  also  reduced.  Mechanical  torque  is 
analogous  to  electrical  current.  Electrical  conductors  including  wires,  integrated  circuits,  motor 
windings,  and  printed  circuit  board  traces  have  electrical  resistance.  That  resistance  restricts  the  flow 
of  electricity  and  causes  a  reduction  in  efficiency  in  the  form  of  its  current  squared  multiplied  by  the 
resistance.  This  energy  is  lost  as  heat  in  the  conductors  themselves.  Careful  analysis  of  the  specific 
conditions  present  in  a  given  system  may  permit  the  use  of  smaller  conductors,  motors,  electronics, 
etc. 
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CONCLUSIONS 

The  benefits  of  implementing  the  presented  algorithm  are  clear.  For  a  minimal  overhead  in 
calculation,  the  user  can  benefit  from  lower  peak  power  use  while  still  reaching  the  given  target  in  the 
least  amount  of  time.  The  time  may  even  decrease  slightly  if  settling  times  are  reduced  by  the 
reduction  in  imparted  forces  caused  by  the  acceleration  decreases.  These  calculations  are  not 
limited  to  laying  of  weapon  systems.  They  are  equally  applicable  to  other  two  axis  systems  such  as 
gantry  tables  and  pan  tilt  devices,  etc. 
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Joshua  Staoo  &  Matthew Tomik _ 

AurhorfPro] «t  Engi  neor  Report  number  (to  he  a$& igned  by  LCSD) 


3061 _  95 _  RDARWSW4 

Extension  Building  Author's/Praject  Engineers  Office 

[Division,  Laboratory,  Symbol) 

PART  1  Must  bo  signed  before  tfte  report  cwi  be  editod- 

a.  The  draft  copy  of  this  report  tips  bean  reviewed  for  technical  accuracy  and  is  approved 
for  editing . 

b.  Use  Distri  button  Statement  A  X  .  B  C_ _ ,  D _ ,  E _ ,  F _ or  X _ for  the 

reason  checked  on  the  continuation  of  thJe  form,  Reason:  _ 

1.  If  Statement  A  Js  selected,  the  report  will  be  reteased  to  the  National  Technical 

information  Service  (NTIS)  for  sale  to  the  general  public.  Only  undassified  r&porte 
whoso  distribution  is  not  limited  or  control  lad  in  any  way  am  released  to  NTIS. 

■  2-  If  Btatemenl  B,  G,  □,  B,  F,  or  X  s  Selected,  the  report  wtll  be  released  to  the  Defense 
Technical  Irtformatton  Center  (DTIC-J  which  will  limit  distribution  according  to  the 
conditions  indicated  in  the  statement 


c.  The  dtetobulten  IJsL  tor  frifc  report  has  b&*n  reviewed  tor  accuracy  and  completeness. 


(  Jose  Santiago 

^IvBiorrtfhtef  ^ 

ytAA'O  lPl<, 

(Data) 

PART  2  To  be  signed  either  when  draft  report  is  submitted  qt  after  review  of  reproduction  copy. 

This  report  is  approved  for  publication. 

John  Gravina 

iSfit/6 

^^ftiwisjon  Chief 

(Date) 

Andrew  Pskowski 

<i '/rtfs' 

RDAR-CJS 

(Oate) 
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